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INTRODUCTION AND RESULTS 

In paper [1] we proposed a general method which allows 
to determine spectra of quantum systems given in the 
Bargmann representation. This representation is very 
useful for systems with a Hilbert space which is a product 
of a finite and infinite-dimensional spaces. A typical 
example is a Hamiltonian for spin degrees of freedom of 
some particles, characterised by n-dimensional a matrices, 
coupled to a bosonic field via annihilation and creation 
operators a, and a^. Then n-component wave function 
ijj = ('01,... ,'0„) is an element of Hilbert space Jf" = 
Jf X • • • X df, where 5f is the Bargmann-Fock Hilbert space 
of entire functions. The scalar product in TC is given by 

{f,9) = - / dxdy, z = x + iy. 

Since operators a and are represented by and multi¬ 
plication by z, respectively, for clearly [dz,z\ = 1, station¬ 
ary Schrodinger equation 11/) = E/ becomes the system 
of linear ordinary differential equations. Usually such 
systems have a certain number of regular singular points 
in the complex z plane, and possibly an irregular point 
at infinity. Determination of the spectrum consists in 
finding such values of E that all the components of wave 
function ipi, with i = 1,..., n, are elements of the Hilbert 
space TC. 

In paper [T] we have shown that a solution of the eigen¬ 
value problem can be reduced to checking the following 
three conditions. 


2. Global conditions. At each singular point we can 
select a local holomorphic solution in such a way 
that they are analytic continuations of each other. 
That is, they are local representations of an entire 
solution. 

3. Normalisation conditions. The entire solution se¬ 
lected above must have a finite Bargmann norm. 

The first two conditions were analysed in detail in [T] 
and in [5] we gave their effective application for determi¬ 
nation of full spectrum of the Rabi system. 

The third condition is highly non-trivial. This is so 
because “most” of entire functions do not belong to fff. 
In order to analyse this condition we need to characterise 
the growth of an entire solution /(z) of a system of linear 
differential equations. In a neighbourhood of a regular 
singular point the growth of /(z) is polynomial. Hence, 
if the considered system has only regular singular points, 
and /(z) is its entire solution, then /(z) has a finite 
Bargmann norm. Normalisation conditions can give non¬ 
trivial constrains on entire solution only if the linear 
system has an irregular singularity, for example at infinity. 

The growth of entire function /(z) is described by 
means of the following function: 

M(r) := max|/(z)|. (1) 

\z\=r 

It is used to define two numbers which characterise prop¬ 
erties of the growth. The order, or the growth order g of 
/(z) is defined as the limit 


1. Local conditions. At each regular singular point 
z = s there exists at least one solution which is 
holomorphic in an open set containing s. 


ln(ln M(r)) 
g := hm sup-;- 

r—>oo Inr 


with 0 < p < oo. 


( 2 ) 


For an entire function of finite order p < oo, its type a is 
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defined as 


the form 


a := lim sup 


In M(r) 


rpQ 


(3) 


If /(z) belongs to 5f, then one can prove the following 
facts 



+ yV'i + 5^2 + 
z'fp2 - yV'2 + g^'i + gzilJi 


Ei)i, 

( 6 ) 

E^2- 


1. f{z) is of order g <2. 

2. li Q = 2, then f{z) is of type cr < f. 

If ^3 = 2 and a = ^, then the question whether /(z) G Jf 
requires a separate investigation. For additional details 

see H]. 

It is well known that in a neighbourhood of a regular 
singular point, a formal procedure allows to find formal 
series which satisfy the equation. It appears that this 
formal series is convergent, so it is a local solution of the 
equation. In a vicinity of irregular singular point a for¬ 
mal procedure gives formal expressions which satisfy the 
equation, however the series are generally divergent. It is 
known that these formal expressions, give the asymptotic 
expansion of a solution valid in in certain sectors radially 
extending from the vertex localised at this singularity. 
Borders of sectors are determined by the so-called Stokes 
lines. The asymptotic expansions change when we pass 
from one sector to another. These changes are governed 
by the so-called Stokes matrices. 

The normalisation condition of an entire solution /(z) 
implies that at each sector S the following integral 


In our paper [T] we calculated the spectrum of H for 
the generic case, i.e., when the above system has two 
regular singular points. This requires that t/^ ^ 4w^, and 
the quantisation of the energy spectrum results from the 
condition that system ([^ admits an entire solution. In 
that case, all entire solutions of have a finite norm. 


In this paper we investigate the remaining cases for 
which = 4tLi^. For U = -|-2a; system can be rewrit¬ 
ten in the form 


■^'liz) 

^ 2 ( 2 ) 


= - zipiiz) H--- 1P2{Z), 

25 

4a;z^ + 2E - ujo 
= -^- 'tpiiz)- 

25 

z -I- uj{2E + Wo)] I ! 


(7) 


The system for the case of 1/ = —2a; can be obtained 
from the above equations by a simple change a;o —> — wq, 
and the interchange ipi with '02- Hence, we consider only 
the case U = +2oj. 


/ |/(z)pe-l^l'dxdg, (4) 

JS 

has a finite value. Thus, at each sector it has good asymp¬ 
totics which guarantees the convergence of the above 
integral. This implies that although a solution /(z) can 
change its asymptotic expansion from sector to sector, 
these changes are restricted only to good asymptotics. 
This property implies that all Stokes matrices have a 
common invariant subspace, in particular case just one 
common eigenvector. 

The aim of this paper is to show that the above consid¬ 
erations can be applied effectively. We consider system 
given by the following Hamiltonian 

H = (^uj+ -I- gaj;{a^ + a), (5) 


System Q has no singular points in the finite part of 
C, so all its solutions are entire functions. Infinity is the 
only and irregular singular point. The system can have 
solutions which grow fast enough to make the Bargmann 
norm infinite. Thus, to determine the spectrum of the 
problem we have to find all values of E for which system 
admits a solution with a moderate growth at infinity, such 
that its Bargmann norm is finite. 

We give a full answer to this question. That is we 
specify explicitly a countable number of energy values 
for which the corresponding eigenfunctions are also given 
explicitly and have a finite Bargmann norm. 

The energy axis is divided into two disjoint intervals, 
each of them with its own countable family of eigenvalues. 
To be more precise, our main result is as follows. Let 


where ax, az are the Pauli spin matrices and wq, oj, g 
and U are parameters. For U = 0, it coincides with the 
Hamiltonian of the Rabi model [S] describing interaction 
of a two-level atom with a single harmonic mode of the 
electromagnetic field. Hamiltonian ([^ was proposed in 
iE]. The term interpreted as a nonlinear 

coupling between the atom and the cavity. 

In Bargmann-Fock representation, the stationary 
Schrodinger equation Eltjj = Etf;, with ip = have 


X 1 + Log ‘^{E + loo/2), ( 8 ) 


be an auxiliary spectral parameter. Then the Hamiltonian 
([^ has entire solutions with a finite Bargmann norm if 
and only if one of the following conditions is fulfilled. 
Either x > I, and 


(x - I) (a;(a; - uJo)+g‘^{x 
uP-^x"^ — 1 


D) 


— 2/1 + 1, 


(9) 
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or a; < — 1, and 

{x - 1) (a;(w - Wo) + g'^{x - 1)) 
Lo'^Vx'^ — 1 


-(2n + l), (10) 


Let us suppose that the system in question can be 
written as 

z~‘^ip'{z) = A{z)ip{z), ipi^z) G C”, (16) 


where n G N. The respective eigenfunctions 


are given by 

= Cexp (-/3±z^) Hn 


( 11 ) 


where Hn{z) denotes the Hermite polynomial of degree 
n, and 


where A(z) is n x n matrix holomorphic in a neighbour¬ 
hood of infinity and Aq := A(oo) is nonzero. Then in 
every sufficiently narrow sector S, the equation has a 
fundamental matrix of solutions of the form 


^{z) = Y{z)z^e^^^\ (17) 

where G is a constant n x n matrix, Q{z) is a diagonal 
matrix whose elements are polynomials of z^^^ with a 
suitable positive integer p, and Y(z) admits an asymptotic 
power series in z~^^p. If we further assume that the 
eigenvalues Xi of Aq, are all distinct, then 


^ x± '/x'^ — 1 

These functions have the growth order g = 2 and type 
cr = /3±. 


ASYMPTOTIC EXPANSIONS 

An asymptotic expansion of a function ijj{z) for z —> oo 
is denoted in the following way 

'ip{z) ~ Oq H-1—^ -I- ... . (12) 

z z^ 

By definition it means that in a given sector 

5(0,5) := {z G C I a < arg(z) < b, |z| > ro}, 
we have 

lim I ^/>(z) - V QjZ-^ ] z'^ =0, (13) 

V 1=0 / 

for an arbitrary m G N, and z G S{a,b). Such series needs 
not have a positive radius of convergence. 

Note that if 

V>(^)e-^(^)~ao + G(z-i), (14) 

for a certain function P{z), then, in particular 

lim '!/;(z)e“^^^^ = Co, (15) 

2—>-00 


• G is a diagonal matrix, 

• p = 1 and Q(z) is also diagonal with leading terms 
of Az«+V(g+1), 

• S can be any sector of central angle not exceeding 
Tr/{q + l). 

The fact that asymptotic expansions are unique allows 
for their formal calculation, and the above tells us, that 
there are always solutions which grow according to the 
prescribed formal expressions. The problem is that for a 
true (not just formal) solution the asymptotic behaviour 
changes from sector to sector—only in rare cases are the 
expansions valid for all values of the argument. 

It can be shown, see [S], that if a function is single 
valued near infinity, its asymptotic expansion is valid for 
all values of arg(z) if and only if it is analytic at z = oo 
and the series is actually convergent . 

In our case the requirement of finite Bargmann norm 
is that 

|^(z)|=G(e"l^l‘'), 0<2, (18) 

and |cr| < 4 when g = 2. So the determination of the 
exponential factor will be crucial, and we shall see that 
the considered model presents us with functions of growth 
order g = 2 and varying types ct, which require detailed 
analysis. 


SPECTRUM DETERMINATION 


and |'0(-z)| = G(exp |P(z)|). 

The general form of solutions around irregular singu¬ 
lar points is known, and there are several classical theo¬ 
rems on them. We choose two, particularly suitable here, 
namely Theorems 12.3 and 19.1 given in [5], which tell as 
the following. 


In the case U = +2oj the considered model is described 
by 0. This system has the standard form ([T^ with 
q = 2. Since the leading matrix Aq is nilpotent, it is 
convenient to perform a shearing transform 


^i{z) = z* ^^^(z), f = 1,2, 


(19) 
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after which the system becomes 



-1 

2uj g‘^(x-l)-ujujo 

g gu:z 2 




^(z), ( 20 ) 


where ^ = [<^ 1 ,^ 2 ]^, a: is a new spectral parameter defined 
by The rank of the system is thus g + 1 = 2, and 
the matrix ^0 has eigenvalues x ± We can 

immediately discard the case of them being equal, for 
then, the system Q is explicitly solvable. For example, 
when cc = 1 we have 

'ipi(z) = ip 2 {z) = - 3wo-z), 

35 

and these functions have infinite Bargmann norm. 

We are thus left only with the simpler case of two 
distinct eigenvalues and we know what form of asymptotic 
expansions to expect. Following Birkhoff [9], let us assume 
the following Laplace integral representation for the vector 
of solutions 

Mz) = ^ exp{z'^y){r]o4y) + ziliAy))‘^y^ (21) 

where i = 1,2. The usual inverse Laplace transform would 
be 


— / exp(zy)? 7 (?/)dy, ( 22 ) 

J c 

However, such 77 would, in general, have essential singu¬ 
larities in the finite complex plane because we are dealing 
with an equation of rank 2. To remedy this, the expo¬ 
nential part needs to be modified in anticipation of the 
growth order 2. This procedure is applicable whenever 
the matrix Aq has different eigenvalues, regardless of the 
rank of the equation, for further details see [9]. 

The contour of integration C is chosen so that the 
integrand assumes the same value at both ends. It could 
be a loop or an open contour extending to infinity. Here, 
the latter is more viable, and C will be one of the curves 
shown in Figure This will allow to integrate by parts 
as follows 


/ exp{z^y)z'^u{y)dy = 

Jc 

[exp(z2j/)z""^u(j/)] c “ y exp( 2 ;^y)z"'"^M'( 2 /)d 2 /, (23) 

and taking the bracket term to be zero. Substituting 
the Laplace representation (21) in equation ([^ we obtain 
a new system for with f = 1 , 2 , by equating 

the integrands to zero. That this is permissible can be 
checked later where the particular rj are recovered but 
we are mostly interested in the asymptotic expansions 
in which case the expressions obtained will be formal 
anyway. 


The integrands are polynomials of z and integrating by 
parts as in (23) we can reduce their degree to at most 1. 
Then the coefficients of z^ and z^ must each be zero so 
that we end up with a system of four equations 


w??ii -f gix - 1)7702 + u}{ 2 y + 1)77(1 = 0 > 
uj{ 2 y + 1)7701 - gix - 1)7712 = 0 , 
w( 2 w 77 oi - g{ 2 y + 2 x - 1)77(2 “ OT12) 

+(wow - g'^ix - l)) 77 oi = 0, 
(woW - g'^ix - l))i]ii + 

w(2w?7(i -f g( 2 y + 2 x - 1)7702) = 0 , 


where prime denotes differentiation with respect to y. 
First two equations can be immediately solved for two 
components 


UJ 


V02 — -7-TT 

5 ( 2 ;- 1 ) 

w( 27/ -f 1 ) 

m 2 - - 7 T- 

gix-i) 


(( 27 /-H 1)77(1-f 7711), 


7 ? 01 - 


(25) 


This leaves two first order equations for the remaining 
components 


( 47 /^ -h dxy + 1 ) 77(0 = {dp\/x‘^ - 1 - 67 / - 3a:) 7710 , 

( 4 //^ -I- 4x7/ -I- 1)77(1 = ^4/7\/— 1 — 27 / — x) 7711, ( 26 ) 


where 

_ (x - 1) (a;(a; - Wo) + ff^(a: - 1)) 
^ 4a;2-\/x^ — 1 


(27) 


Although the particular equations decouple, this is obvi¬ 
ously not a generic feature. It makes sense then to treat 
them as a system with two regular singular points at 

ai ,2 = ^ (-a: ± Va;^ - l) , (28) 


and one at infinity. The characteristic exponents are 


h & {p - \ , P - \} at oi, 

P2& {-p-l,-p-\] at 02, 

/?oo G {§,\) at 00 . 


(29) 


In general, the behaviour of a solution as it is continued 
around all the singular points cannot be specified explicitly 
as the monodromy group of the equation is unknown. 
Fortunately, solutions of system ( [^ are given explicitly 

? 7 oi =biiy- aif-^/'^iy - 02)"'’"^^^, 

7711 =62(y-al)^-l/"(7/-02)-'’-'/^ 

where 61^2 are arbitrary constants. The explicit forms of 
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these solutions will give us full information about which 
contours to choose and how the continuation of any solu¬ 
tion behaves. 

Each singular point of the solution r]i j{y) corresponds 
to a different exponential factor in the asymptotic ex¬ 
pansion, which can be seen as follows. Let us look at 
the generic case first, when the contour can be chosen 
as path Ci that runs from infinity to the vicinity of the 
point y = Oi, around it in the positive direction, then 
back to infinity, as depicted in Figure This requires 





FIG. 1. Contours. ai = 5 — l/(4a2) 


that the characteristic exponent at that singularity at is 
not an integer, for otherwise the integral vanishes. Since, 


for the system (26) for all singular points the differences 


of exponents is 1/2, one solution is always suitable. 


The Laplace integral (21) for any function ri{y) which 


has a local power series around y = a with characteristic 
exponent 13 may be rewritten as 


27ri 


'C 


°‘\y - aY^aniy - a)'^Ay 


n=0 

OO 


V(-R - ■ 


^2(/3+i) r(—/3 — n)z^'^ ’ 

n=0 


(31) 


will be dealt with separately, below. 

It is important to keep in mind that the contour had 
to be chosen so that the integrand vanishes as y —> oo 
and that the integral (21) is convergent - this is ensured 
by Ke{z^y) < 0. In other words, for a given argument of 
z there is just a sector of the complex plane of y where 
the integral converges and vice versa: for a given ray 
(arg(j/ —a)) in the y plane, there is a region of convergence 
in the z plane. An example is shown in FigureBecause 
the components y grow at most with some power of y at 
infinity, the integral (21) converges fast enough, for the 
function to be holomorphic and it is indeed permissible 
to differentiate it as we have done. 

To extend the function outside one sector, the integral 
can be analytically continued by deforming the contour, 
rotating it around a given singularity. Because the contour 
has to pass through the other singularity, the solution 
associated with one asymptotics could transform into 
a combination of both, as shown in figure This is 
the Stokes phenomenon, and the argument of ai — 02 
determines the so-called Stokes lines in the z plane, which 
is customarily given by the condition Re(z^(ai — 02 )) = 0; 
crossing the line causes the asymptotic expansions to mix. 
There is the freedom of choosing the starting contour, 
which corresponds to choosing a particular solution, and 
the first Stokes line where the deformation is required 
depends on that choice. But if we continue a full basis 
increasing the argument of z we eventually need to take 
into account all of the lines. For the considered case, 
ai — 02 is real, so there are two Stokes lines which pass 
through the origin and form angles ±7r/4 with the x-axis. 


Upon crossing, the original integral decomposes into in¬ 
tegrals over new contours, which could be written shortly 
as 


A [ = [ + [ - [ , (32) 

JCi JC- JC[ JC- 

where A stands for analytic continuation, and primes are 
used to distinguish between contours in the two regions: 
Im(?/) > 0 and Im(y) < 0, although their shape remains 
the same. This could be written using the Stokes matrix 
as 


where we have used the integral representation of the 
gamma function, and omitted indices for brevity. The 
series is not, in general, convergent in the whole complex 
plane of y, so the integration term by term using the 
gamma function leads to a formal, divergent in most 
cases, asymptotic expansion. It can be shown that this is 
the correct asymptotic series nni, so that it is immediately 
evident that each singular point ai gives rise to one type 
of exponential factor exp(Q;iZ^) in the asymptotic series. 
Note that 40102 = 1 and in view of the previous section 
we require |oi| < 1/2, so only one of the contours will lead 
to the desired growth type. The special case of |oi| = 1/2 


A{ei,eY = (e', e') , (33) 

where is the solution obtained with Ci, hence with 
fixed asymptotics of exp(aiZ^), and likewise for e' and 
C[. By convention, the matrix on the right-hand side is 
the inverse of the Stokes matrix. 

As (3i is not an integer, whether the additional integral 
actually contributes, depends on the characteristic expo¬ 
nent of the solution around aj. If /3j is not an integer 
either, the Stokes multiplier is obviously not zero. On the 
other hand, for an integer /3j the contour Cj can simply 
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ImCy) 


Im(z) 


Re(y) 


Re(z) 




FIG. 2. The corresponding sectors in the y and 2 ; planes for 20 = 1 + i/2. Note that both shaded regions on the right fulfil 
Re(y 2 ^) < 0 for any given y with arg( 2 /) = tp. 


Im(t/) 



Re(j/) 


FIG. 3. Deformation. 


then 


A{ei,ej) = (e',e') . (35) 

For completeness we note that the contour Lj can be 
deformed analogously to the above, this time winding 
around ai. Because Pi is not an integer, the contour 
cannot pass through the singular point and the additional 
path will be exactly C^. In terms of the Stokes matrix we 
have 


^(e',e') = (e",e")(j j) • (36) 

In other words, the asymptotics of the solution re¬ 
mains unchanged, while the asymptotic series ej actually 
exhibits the Stokes phenomenon because Pi is not an 
integer. The double primes signify contours identical to 
the original ones for Im( 2 /) > 0 but the solutions are not 
single-valued in general so e" are some multiples of e^. 


be deformed through the singularity so the integral is 
identically zero, because 



This does not mean that the asymptotic solution van¬ 
ishes, but simply that the contour Lj , which is a half line 
from aj to infinity, has to be used instead of Cj to define 
the other asymptotic solution ej. The Stokes matrix is 


Note, that when Pj is an integer. Pi cannot be, since 
their sum is —3/2 or —1/2 for the two solutions of (30), 
respectively. Thus, the conditions that the contour is 
suited for the noninteger exponent at ai and that the 
Stokes multiplier vanishes are consistent. There is no 
need then to investigate the case of integer Pi, because for 
any p we always have at least one solution with noninteger 
exponent at each of the singular points. It follows that 
we always have two different solutions with different 
asymptotic expansions obtained by means of the contours 
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C and, as stated in the previous section, the fundamental 
solutions have to be expressible by means of those two 
Ci in any sector. In each case then, we know how both 
asymptotic solutions mix when crossing the Stokes lines 
and hence, how the growth of any actual solution changes. 
For the system under consideration the situation is very 
simple, because only the asymptotic series with smaller 
Qfi is suitable, so that all of its Stokes multipliers must be 
zero to ensure finite norm. 

It remains to verify that when x'^ < 1 and hence |(ai| = 
I /2 the resulting solutions are not normalisable, in close 
analogy to the x = 1 case, where the components of the 
wave function behave like The exponential factor 

in the present case can be written as 

exp , for (37) 


and since the equation determining aj is real we have 
Oil = 0(1 so 6*2 = — 6*1 ^ 0. Together with the Bargmann 
measure this results in the leading factor of the form 



rdrdip, 


(38) 


where we took z = re"^. The integral over r converges pro¬ 
vided that the sine does not vanish, i.e. in sectors that do 
not include the directions (p + 0jl2 = kir, k G Z. If there is 
no Stokes phenomenon for the solutions with asymptotic 
factors of e“’^ , it follows that the norm is not finite. If 
the two solutions could interchange, we would necessarily 
have to construct a solution which has the asymptotics 
of exp(aiz^) in the vicinity oi ip + 6 * 2/2 = 0. But then, 
continuing it through the angle of tt in arg( 2 ;) = ip, we 
would end up in a sector containing the next problematic 
line (/? -|- ^ 2/2 = TT. The corresponding contour in the y 
plane would have to rotate through 27r and the Stokes 
phenomenon would cause the asymptotic expansion to 
now contain also the series with exp{a 2 Z^) whose norm 
would diverge in the second sector. Thus regardless of 
the Stokes multipliers, although the solutions are both 
entire, none of them is normalisable for — I < a; < I. One 
could interpret this as a continuum spectrum with gener¬ 
alised eigenstates, because to each entire function in the 
Bargmann space corresponds a well behaved distribution 
solution of the original problem in the (K) space [3] . 

Taking the above into account, we can finally give 
conditions for the original entire function to be normal¬ 
isable. First, for the desired asymptotics one must have 
\ai\ < 1/2 so that the solution starts as normalisable 
within some Stokes sector, and second f]j G N, so that 
it remains such, as it is continued throughout the whole 
complex plane. The second condition is going to be of 
the form ±p — 3/4 € N or ±p — 1/4 S N as there are two 
solutions of ( [^ to choose from, but these can be written 
together as ±4p = 2n -I-1, n G N. 


Returning to the spectral parameter x, the normalisa¬ 
tion, and quantisation, conditions read: 


|ai| < ^ a; > 1 , 

—4p = 2n+ -1-1, n_|_ G N, 


laal < ^ X < -1, 

4p = 2n_ -1-1, n_ G N. 


(39) 

(40) 


or, equivalently. 


sgn(a:) 


{x - 1) {uj{u} - Wo) -I- g‘^{x - 1)) 
oj'^y/x'^ — 1 


= 2n + l, (41) 


for > 1 and n G N. The energies of these spectra 
are presented in Fig. |^and[^ note that the lower levels 
become denser as E gets closer to the continuous (non- 
normalisable) region a;^ < 1 . 

It is evident, that with these conditions there will always 
be one characteristic exponent that is a negative integer. 
In this case one of the contours C will be equivalent to 
a circle and the corresponding Laplace integral can be 
immediately evaluated using Cauchy’s theorem. It will 
yield a finite number of terms so that the solution will be 
a polynomial in z with an exponential factor. A change 
of variables then reveals that each interval of x has its 
own sets of eigenstates expressible in terms of Hermite 
polynomials, namely 


li^iiz) = Cl exp 
when X < — 1 , or 
= C2exp 


—x — a/x^ — I 


2 " Vx^-lz 


=)' 

(42) 


—X -I- Vx^ — I 






X^ — I Z 


when X > 1 , and in each case the other component is 

Mz) = "" 


gix- I) 


(43) 

S 

(zi/jiiz) + ii^iz)). (44) 


CONNECTION WITH THE WHITTAKER 
EQUATION 


Let us rewrite system § as one equation of the second 
order. It reads 

w" + 2xzw' + + X — my/x^ — I^ w = 0, (45) 












xy 


(47) 



FIG. 4. The upper spectrum for uj = 1 and cjq = 1/2. 
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FIG. 5. The lower spectrum for cu = 1 and loq = 1/2. Thin 
lines denote constant energy levels. 


where w{z) = x is defined by ([^, and 

(x - 1 ) [u{uj - uJo) +g'^ix- 1)1 
m := - 

uP-\/x'^ — 1 

Now, we introduce new independent variable 
change the depended variable putting y = 


(46) 

y, and 
•s/x^ — 1 


and 


x{y) = y^/'^exp 


2^/x^ - 1. 


w{z). 


We obtain the Whittaker equation 

v'\y)-(\-^ + v{y) = 0, (48) 

\4 y y^ J 

with parameters 

K = -^m, f^=\- ( 49 ) 

According to Proposition 2.5 in m, there exists a basis 
of solutions of this equation such that two Stokes matrices 
have the form 


Si 


1 a 
0 1 


and S 2 


1 0 
P 1 ’ 


(50) 


where a and /3 are some complex numbers. Moreover, the 
Stokes multiplier a vanishes if and only if either k + fj,, 
or K — y is a non-negative half integer. Similarly, the 
Stokes multiplier /? vanishes if and only if either —k — /i, 
or — K -|- /r is a non-negative half integer. The reader can 
check, that in the considered case these conditions exactly 
coincide with those given by Q and 

Now, it is interesting to observe that for the Whittaker 
equation if one of the Stokes multipliers vanishes, then it 
is solvable in the sense of differential Galois theory—all 
its solutions are Liouvillian. More precisely, the identity 
component of its differential Galois group is Abelian, see 
Corollary 2.1 in m- We have the same conclusions for 
equation (451 describing our problem. Thus, if E is an 
eigenvalue then; 


1. A Stokes multiplier vanishes; 


2. All solution are Liouvillian; 


3. The identity component of its differential Galois 
group is solvable. 

These observations are very important for application. 
If we are able to prove inverse implications, then we 
have a strong tool for quantisation and finding exactly 
solvable systems. In fact, for low order equations there 
exist effective algorithms for studying their differential 
Galois groups and Liouvillian solutions. 

Concerning the original equation, we have the full char¬ 
acterisation of its local Galois group at infinity. Although 
in this case it coincides with the full global group, the 
conclusions of [12] suggest, that for systems with other 
singular points, the Galois group at infinity alone gives 
quantisation conditions. Namely, for energy values which 
belong to the spectrum, the group of the corresponding 
time-independent Schrddinger equation is solvable. 
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CONCLUSIONS 

[ 2 ] 

According to our knowledge, this is the first example 
of application of conditions deduced from an analysis j-gj 
of the Stokes phenomenon for a quantum system in the 
Bargmann representation. We only know about appli¬ 
cation of the finite norm condition for square-integrable [4] 
functions to determine spectra for one-dimensional quan¬ 
tum systems in position representations. Finiteness of 
this norm implies some conditions on local solutions con- jg, 
taining real ±oo, see e.g. quantisation of quartic oscillator 
in m, and related problem of Stokes phenomena for pro¬ 
late spheroidal wave equation considered in [T^]. In our [7] 
considerations, asymptotic solutions at all sectors, not 
just real infinity, are involved. It appears that in some 
cases these conditions are calculable analytically and then 
one can obtain the spectrum explicitly. 
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